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ABSTRACT 



We report the first sub-arc second (0.65" x 0.51") image of the dimethyl 
ether molecule, (CH 3 ) 2 0, toward the Orion Kleinmann-Low nebula (Orion-KL). 
The observations were carried at 43.4 GHz with the Expanded Very Large Ar- 
ray (EVLA). The distribution of the lower energy transition 61,5 — 60,6, EE (E u 
= 21 K) mapped in this study is in excellent agreement with the published 
dimethyl ether emission maps imaged with a lower resolution. The main emis- 
sion peaks are observed toward the Compact Ridge and Hot Core southwest 
components, at the northern parts of the Compact Ridge and in an intermediate 
position between the Compact Ridge and the Hot Core. A notable result is that 
the distribution of dimethyl ether is very similar to that of another important 
larger O-bearing species, the methyl formate (HCOOCH3), imaged at lower res- 
olution. Our study shows that higher spectral resolution (WIDAR correlator) 
and increased spectral coverage provided by the EVLA offer new possibilities for 
imaging complex molecular species. The sensitivity improvement and the other 
EVLA improvements make this instrument well suited for high sensitivity, high 
angular resolution, molecular line imaging. 

Subject headings: ISM: individual objects — ISM: jets and outflows — ISM: 
molecules 



1. Introduction 

The Orion complex is among the most studied sources in our Galaxy. At a distance of 



420 pc (see, e.g. Sandstrom et al. 2007 Hirota et al. 2007 Menten et al. 2007) this is the 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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nearest site of recent high mass star formation. The prominent H II regions NGC 1976 or M42 
and NGC 1977 are on the near side of a large molecular cloud associated with Lynds dark 



cloud L1640 (see O'Dell 2001). At the rear of the H II region is the Orion Molecular Cloud, 
OMC-1. The brightest molecular emission in OMC-1 is near the H II region NGC 1976 (see, 
e.g. Tatematsu et al.||l993 ). Plume et al. (2000) measured the J=5-4 transition of 13 CO over 
a 0.5° by 2° area with a 3.2' resolution. The brightest emission has an extent of ~15' in the 
north-south by ~5' in the east-west direction. Within this region of extended emission is a 



prominent maximum, the Kleinmann-Low (KL) infrared nebula (see, e.g. Dougados et al. 



1993 Gezari et al. 1998) which exhibits emission from complex molecules, masers, outflows 
and warm dust. In addition to extended IR emission from warm dust, there are ~20 compact 



near-IR sources found toward Orion-KL (Dougados et al. 1993) 



In quasi-thermal continuum emission from dust grains, Orion-KL shows an overall ex- 
tent of 10" in a by 15" in 5. With higher angular resolution, the form has a "V" shape, 
with the symmetry axis at position angle ~20° east of north. Within this region, there is 
fine structure, with molecular species peaking at different positions. This reflects abundance 
differences, not excitation effects. 



In the methyl formate (HCOOCH3) image (see Fig. 4 of Favre et al. 2011) and in the 



frequency-integrated maps of several molecular tracers (see Fig. 5 of Guelin et al. 2008) 



the molecular distribution is seen as a V-shape with each arm having an observed length 
of 12" by 3", with the major axis at a position angle ~20° east of north having an opening 



angle of 38°. From Plateau de Bure Interferometer (PdBI, Favre et al. 2011) and EVLA 
data (A. Remijan, unpublished), the methyl formate emission taken with resolutions of 
1.8" x 0.8" and ~5" peaks at a J20 oo = 05 h 35 m 14"09, S J2000 = -05°22'36"7. This is the center 
position of the compact ridge region. The compact ridge shows a larger abundance of oxygen- 
containing species. In CH 3 CH 2 CN, ethyl cyanide, the center of the Hot Core peaks at «j2ooo 
= 05 h 35 m 14"6, 5j2ooo = -05°22'29", with observed sizes of 9" by 4" at PA. ~20° east of 



north. This region shows an excess of nitrogen bearing species (see, e.g. Blake et al. 1996 



Friedel fc Snyder||2008[ |Wilson et aL||2000| |Wright et aL|l996[ ). The chemistry of Orion-KL 
is rich, with rather large abundances of complex molecules. There are two possibilities: (1) 
gas-phase reactions between ionized and neutral species, and (2) formation on dust grains, 
followed by liberation from grains by radiation from nearby IR sources. 

In addition to extended dust and molecular emission, Orion-KL also has a number of 



radio continuum sources (Garay et al. 1987 Churchwell et al. 1987 Menten & Reid 1995 



Rodriguez et al. 2005). Accurate high angular resolution radio data show that the highly 



obscured source "I", and the Becklin-Neugebauer object (BN Becklin & Neugebauer 1967) 



show large proper motions (Plambeck et al. 1995). Subsequently, (Rodriguez et al. 2005 
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Gomez et al. 2005, 2008; Goddi et al. 2011) have found proper motions of source I, BN and 



the radio counterpart of source "n" . 

Given the large amount of activity in the Orion-KL region, high angular resolution is 
needed to separate the various influences. There are two coupled problems to be solved. 
The first is the chemistry of the Hot Core and the Compact Ridge. Although the centers 
of these regions are separated by only ~10" (0.02 pc), their chemistry is quite different. 
In addition, models must accurately predict the abundances of complex species and their 
survival in the presence of the activity as shown by the large proper motions of sources BN, 
I and n. Given these questions, the first task is to image complex molecular species with 
arc second or better angular resolutions, to test chemistry models. For this purpose, we 
have imaged the KL nebula in dimethyl ether (CH 3 ) 2 0, which had been mapped with lower 



spatial resolution by Friedel & Snyder (2008) and Guelin et al. (2008). 



In Sect. [2] we present our observations and we briefly describe the data reduction 
methods. Our maps are presented in Sect. [3j We compare our results to previous studies in 
Section |4} Conclusions are presented in the last Section. 



2. Observations 

A quartet of dimethyl ether ((CH 3 ) 2 0) lines (see Table [T| were observed on 9 Nov 2010 



(1 hour session) and 19 Jan 2011 (4 hours session) with the NRAO EVLA (see Perley et al. 



2011) in the "C" configuration (proposal code: 10B-223). Dual polarization observations 
were made toward a single pointing (aj2ooo — 05 h 35 m 14"2, <5j 2 ooo — -05°22'36"0), tuned to 
the frequency of 43433.5449 GHz. Here we discuss the 19 Jan 2011 observations. The 
targeted lines were the 6^5 — 60,6 lines lying near 43.44757 GHz. AE, EA, AA and EE lines 
appear close together in a distinctive "pitchfork" pattern. The upper energy levels are only 
21 K above ground and easily excited. 

Winds were low and a few scattered clouds decorated the sky under dry conditions on 
19 Jan. Using the measured weather parameters and an atmospheric model generated from 
the ATM program within CASA we estimate that the atmospheric opacity was 0.058 during 
the observations; we have corrected the data for that assumed value. 

The full width at half power (FWHP) at the observing frequency of the primary beam 
of each antenna is The WIDAR spectrometer provided 256 channels over a 16 MHz 
bandwidth and a channel spacing of 62.5 kHz (0.43 km s -1 ) in each of the two orthogonal 
polarizations. Centered near the nominal velocity of the ambient cloud, vlsr=8 km s -1 , 
111.2 km s -1 of velocity coverage was available in total. 
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J0541-0541 (measured flux density S 43 .43 = 0.386±0.001 Jy) served as the gain calibra- 
tor for Orion-KL using fast-switching (2 m on source, l m on the calibrator). Data calibra- 
tion followed standard high-frequency procedures in CASA, including using a model for the 
brightness distribution of the absolute flux calibrator J0137+331 (3C48), which gives a flux 
density of 843.43= 0.528 Jy. J0319+4130 (3C84) was used as bandpass calibrator (843.43=13.6 
±0.5 Jy). 

The continuum emission was subtracted from the January dataset and imaged. The 
final continuum and continuum-subtracted images for analysis have been corrected for the 
response of the primary beam. 

The final naturally weighted maps reached a (5a) sensitivity of 20 mJy per beam (39 K) 
in the line images and 3 mJy per beam in the continuum images, as expected. The beamsize 
is 0"65 x 0"51 at position angle 38°. 



3. Results 

3.1. Continuum emission 

The 43 GHz continuum emission observed toward Orion BN/KL is shown in Fig. [T] 
Two continuum sources, the BN object and the radio source I, are detected above 5cr with 
a peak flux density per synthesized beam S u , of 22.3 mJy beam -1 and 9.8 mJy beam -1 , 
respectively. Some of the structure of the Hot Core source is resolved. Emission from the 
BN object also appears quite compact. The total flux densities S^, are given in Table [2j 



3.2. Dimethyl ether (CH 3 ) 2 

The dimethyl ether emission map at 43.447 GHz is shown in Fig. [2] The observed four 
main molecular emission peaks, hereafter DME1, DME2, DME3 and DME4, are located 
toward the Compact Ridge, the Hot Core southwest (Hot Core-SW), in an intermediate 
position between the Compact Ridge and the Hot Core-SW and in the north of the Compact 
Ridge, respectively. Most of the peaks appear with a LSR velocity of 7.4 km s -1 . As for 
many O-bearing species the emission distribution presents an extended V-shape linking the 



radio source I to the BN object (see, e.g. Guelin et al. 2008) 



The high spectral resolution provided by the WIDAR correlator allows us to distinguish 
between the AA and EE transitions, although the AE and EA transitions are blended (see 
Table [TJ. Spectra of the lower energy transitions 6^5 — 6o,6, AA, EE and AE/EA (E u = 
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21 K) observed toward the regions DME1 to DME4 are presented in Fig. [2j The observed 
line parameters (velocity, intensity) of the detected (CH 3 ) 2 6^5 — 6 ,6, EE transition are 
summarized in Table [U 



4. Comparison with other results and analysis 
4.1. Comparison with continuum maps 



The total flux densities S u (see Table |2J), obtained toward source I, in the Hot Core 
region, and the BN object are in very good agreement with the 43 GHz Very Large Array 



(VLA) values reported by Menten fc Reid (1995); Chandler fc Wood (1997); Reid et al. 



(2007); Goddi et al. (2011) 



4.2. Single dish observations 

Green Bank Telescope (GBT) molecular surveys have been undertaken from 42.3 GHz 
to 43.6 GHz ( |Goddi et al-pOll] ) and from 42.7 GHz to 45.6 GHz (Wootten, unpublished and 
Guelin et al. 2008). A comparison between these surveys and our EVLA spectra has allowed 
us to conclude to within the uncertainties, that all flux measured in either GBT survey for 
the (CH 3 ) 2 0, 61,5 — 60,6 transition is recovered in our EVLA images. 

The excitation of the dimethyl ether emission may be estimated from fitting a rotational 



temperature to the three lines spanning an energy range 21-147 K from Goddi et al. (2009) 



or from two lines spanning the same range from Guelin et al. (2008). From these data, line 



sets observed simultaneously on the GBT, we estimate T rot ((CH3) 2 O)~120 K, in agreement 
with estimates from transitions of HCOOCH3 of Favre et al. ( |2011 ). Note that the brightness 
temperature of the dimethyl ether line reported here toward the peak of the compact ridge 
(DME1, see Fig. |2| is ~140 K and that this line is optically thick (r estimated to be 2.9), 
consistent with the rotational temperature in the broader (16") GBT beam. Using lines 



of acetone (CH 3 ) 2 CO in the same GBT spectrum, Goddi et al.| (2009) found a rotational 
temperature twice this value, suggesting an origin in the Hot Core. Higher resolution images 
are needed to understand how the distributions of otherwise similar O-bearing molecules 
differ, and to glean lessons on how large molecules may form and vanish in molecular clouds. 
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4.3. Interferometric observations 



Using the Combined Array for Research in Milimeter-Wave Astronomy (CARMA), 



Friedel h Snyder (2008) have imaged the dimethyl ether transition 13 ,i3 — 12i,i2 (E u 



86 K) with a beam size of 2.5" x 0.85". Our observed (CH 3 ) 2 emission peaks (DME1, 
DME2 and DME4, see Fig. [2]) are present in their lower resolution data Q CARMA obser- 
vations reveal a similar LSR velocity at 7.6 km s _1 . However, their lower spectral resolution 
(by a factor of 1.4 compared to that of ours) did not allow them to resolve the AA, EE and 
AE/EA lines. 

A comparison between our EVLA map and the Plateau de Bure Interferometer (PdBI) 
dimethyl ether transition 122,10 — 113,9 (E u = 78 K) imag ed by |Guelin et al7] fl2008| fj with 
a synthesized beam of 3.3" x 1.7", shows very good agreement. With a higher- angular 
resolution the distribution emission of the (CH 3 ) 2 mapped with EVLA presents the same 
main strongest emission peaks within the entire region: DME1, DME2 and DME4 (regions 



A, B and C, respectively in Fig. 5 of Guelin et al. 2008). 



Our EVLA sub-arc second (0.65" x 0.51") image of the dimethyl ether clearly reveals a 
new emission peak DME3, located between the Hot Core-SW and the Compact Ridge (see 
Fig. |2l. 



4.4. Relation to methyl formate (HCOOCH 3 ) and shocks 



There is a debate whether dimethyl ether is formed on grain surfaces or in the gas phase 



via a path involving methanol as a precursor (Garrod et al. 2008 Peeters et al. 2006). This 



largest O-bearing molecule is highly abundant in star formation region, particularly in hot 
cores ( Nummelin et al.||2000 Sutton et al.||1995" Turner|[l9 91). Studies of the astrochemistry 
of this molecule, as well as its spatial distribution will bring important data for understanding 
high-mass star forming regions. 

A comparison between the dimethyl ether distribution with that of another important 
O-bearing molecule, methyl formate HCOOCH3 (see, e.g. Favre et al.||2011 ) shows that both 
of these species peak at the same locations within the same velocity range. In particular, 



1 Note that the Compact Ridge position defined in 



Friedel & Snyder 



used in this study (e.g. reference taken from Beuther et al. 2005). Their IRc5 position is closer to our 
Compact Ridge position (1.5" away). 



(2008) is not the the same as that 



2 Note that a misprint with the dimethyl ether quantum numbers appears in the Table 1 of 



Guelin et al. 



(2008). 
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these share common emission peaks in the north of the Compact Ridge, in the Compact 
Ridge itself, toward the Hot Core-SW and toward an intermediate region linking the Hot 
Core to the Compact Ridge. These four emission regions are shown in Fig. [2] (see labels 
DME1 to DME4 in the present letter and MF1 to MF4/5 in |Favre et aL]|2011[ ). Hence the 
formation of these two molecules must have some relation, regardless whether the mechanism 
involves gas phase or grain surface formation. 



Favre et al. (2011) have shown that shocks could be responsible for the HCOOCH3 pro- 



duction in the Compact Ridge region. In particular, the clear association observed between 



the emissions of the 2.12 /mi vibrationally excited H2 (Lacombe et al. 2004) and the methyl 



formate could result from shocks between the interstellar material and bullets owing to the 



close encounter of the sources I, n and BN 500 years ago (Zapata et al. 2009; Bally et al. 



2011). The similar spatial distribution observed suggests that the same mechanisms could 



be at the origin of the release of (CHa^O (itself or a precursor) from ice grain mantles. 



5. Conclusion 

We have imaged the distribution of the dimethyl ether 6i j5 — 6 0i 6, EE (E u = 21 K) tran- 
sition with the high-angular resolution of 0.65" x 0.51" using the EVLA. The high resolution 
provided by the WIDAR correlator, allowed us to spectrally separate transitions AA, EE 
and AE/EA. Our study shows the great capabilities offered by this new interferometer at 
such high spatial and spectral resolution. 

Our results confirm the three main emission peaks of (CH 3 ) 2 observed in the previous 
studies with a lower angular resolution. A new dimethyl ether emission peak has been 
identified thanks to the high spatial resolution of these EVLA observations. The similarity 
between spatial distributions of this molecule and HCOOCH3 and vibrationally excited H 2 
toward the Compact Ridge suggest that same production mechanisms that is, shocks may 
be involved. 

We are grateful to the entire EVLA staff who produced such an excellent instrument. 

Facilities: The Expanded Very Large Array of the National Radio Astronomy Observa- 
tory. 
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Fig. 1. — 43 GHz continuum emission observed with the EVLA toward Orion BN/KL. The 
first contour and the level step are 4.7 mJy beam -1 and the clean beam (0.65" x 0.51") is 
shown in the bottom left corner. Position of the BN object is (aj 2 ooo = 05 h 35 m 14!1094, 5j2ooo 
= -05°22'227724) and the position of the radio source I is (a J20 oo = 05 h 35 m 14!5141, 5 J2000 = 
-05°22'307575) QGoddi et aLpOlI). 



13 




Fig. 2. — Dimethyl ether (CHs^O emission map observed with the EVLA toward Orion 
BN/KL. The synthesized beam, 0.65" x 0.51", is shown in red in the lower left corner. 
The red contours show the 43 GHz continuum emission of the BN object and the radio 
source I (see Fig. [I]). Spectra displayed in this figure are observed toward the four main 
(CH 3 ) 2 emission peaks (DME1 to DME4), i.e. in the following directions: the Compact 
Ridge (bottom right), the Hot Core-SW (top left), the intermediate region linking the Hot 
Core to the Compact Ridge (bottom left) and the north of the Compact Ridge (top right). 
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Table 1: Spectroscopic line parameters 3, of dimethyl ether (CH3)20. 



Frequency 13 


Transition 


( Sij/i 2 ) 


E u 


(MHz) 








(D 2 ) 


(K) 


43446.4708(12) 


61,5 


- 6 ,6 


AE 


14.0 


21.1 


43446.4713(12) 


61,5 


- 6 ,6 


EA 


28.1 


21.1 


43447.5415(11) 


61,5 


- 6 ,6 


EE 


112.2 


21.1 


43448.6120(14) 


6l,5 


- 6 ,6 


AA 


42.1 


21.1 



"All spectroscopic data from (CHa^O taken from Groner et al. (1998) available from the JPL molecular line 
catalog ( Pickett et al.(l998 ) at Splatalogue (www.splatalogue.net, Remijan et al.|[2007 ). 
b Errors are 2a. 



Table 2: Total flux densities obtained at 43 GHz for the sources BN and I. 



Source 


«J2000 a 


^J2000 a 


Uncertainty 13 






m 


(° ' '0 


(mJy) 


(mJy) 


BN 


05 35 14.1094 


-05 22 22.724 


3 


29 


I 


05 35 14.5141 


-05 22 30.575 


3 


16 



Goddi et al. (2011) 



6 5ct. 



Table 3: Line parameters for the detected (CHa^O 61,5 
toward Orion BN/KL. 



60,6, EE transition at 43.4 GHz 



Region 


Peak label 


V LSR 


Av LS r 


Io 






( km s _1 ) 


( km s -1 ) 


(Jy beam -1 ) 


Compact Ridge (CR) 


DME1 


7.4 


0.9 


64 


Hot Core-SW (HC-SW) 


DME2 


7.4 


1.9 


14 


Intermediate of CR and HC-SW 


DME3 


7.4 


1.7 


35 


Compact Ridge North 


DME4 


6.9 


1.5 


40 



